1.1 Installing R

Go on this link to download R: https://cran.rstudio.com/

Select the version that works for your operating system, and download the latest release (R-3.6.0).

Figure 1.1: Download R.

Once you’ve downloaded R, install it following the instructions on the screen.

1.2 Installing R Studio

Go on this link to download R Studio: https://www.rstudio.com/products/rstudio/download/#download

And then download the version that works for your operating system.

Figure 1.2: Download R Studio.

Once you’ve downloaded R Studio, install it following the instructions on the screen.

3.1 R Studio

• A great integrated development environment (IDE) in which you can do all your R coding.
• let’s change some settings first
• R Studio cheatsheet

Figure 3.1: General preferences.

Make sure that:

• Restore .RData into workspace at startup is unselected
• Save workspace to .RData on exit is set to Never

Figure 3.2: Code window preferences.

Make sure that:

• Soft-wrap R source files is selected

This way you don’t have to scroll horizontally. At the same time, avoid writing long single lines of code. For example, instead of writing code like so:

ggplot(data = diamonds, aes(x = cut, y = price)) +
  stat_summary(fun.y = "mean", geom = "bar", color = "black", fill = "lightblue", width = 0.85) +
  stat_summary(fun.data = "mean_cl_boot", geom = "linerange", size = 1.5) +
  labs(title = "Price as a function of quality of cut", subtitle = "Note: The price is in US dollars", tag = "A", x = "Quality of the cut", y = "Price")

You may want to write it this way instead:

ggplot(data = diamonds, aes(x = cut, y = price)) +
  # display the means
  stat_summary(fun.y = "mean",
               geom = "bar",
               color = "black",
               fill = "lightblue",
               width = 0.85) +
  # display the error bars
  stat_summary(fun.data = "mean_cl_boot",
               geom = "linerange",
               size = 1.5) +
  # change labels
  labs(title = "Price as a function of quality of cut",
       subtitle = "Note: The price is in US dollars", # we might want to change this later
       tag = "A",
       x = "Quality of the cut",
       y = "Price")

This makes it much easier to see what’s going on, and you can easily add comments to individual lines of code.

RStudio makes it easy to write nice code. It figures out where to put the next line of code when you press ENTER. And if things ever get messy, just select the code of interest and hit cmd + i to re-indent the code.

Here are some more tips on how to write nice code in R:

• Advanced R style guide

3.2 Getting help

There are three simple ways to get help in R. You can either put a ? in front of the function you’d like to learn more about, or use the help() function.

?print
help("print")

Tip: To see the help file, hover over a function (or dataset) with the mouse (or select the text) and then press F1.

I recommend using F1 to get to help files – it’s the fastest way!

R help files can sometimes look a little cryptic. Most R help files have the following sections (copied from here):

Title: A one-sentence overview of the function.

Description: An introduction to the high-level objectives of the function, typically about one paragraph long.

Usage: A description of the syntax of the function (in other words, how the function is called). This is where you find all the arguments that you can supply to the function, as well as any default values of these arguments.

Arguments: A description of each argument. Usually this includes a specification of the class (for example, character, numeric, list, and so on). This section is an important one to understand, because arguments are frequently a cause of errors in R.

Details: Extended details about how the function works, provides longer descriptions of the various ways to call the function (if applicable), and a longer discussion of the arguments.

Value: A description of the class of the value returned by the function.

See also: Links to other relevant functions. In most of the R editors, you can click these links to read the Help files for these functions.

Examples: Worked examples of real R code that you can paste into your console and run.

Here is the help file for the print() function:

Figure 3.3: Help file for the print() function.

3.3 Installing and maintaining packages

What makes R powerful is the large number of packages that have been written for R. You can install a new package like so:

install.packages("tidyverse")

You can also install multiple packages at the same time, by concatenating the package names using the c() function:

install.pacakges(c("tidyverse","broom"))

To make sure that your packages remain up to date, you can go to Tools > Check for Package Updates ... in R Studio.

Figure 3.4: Help file for the print() function.

You can then click Select All and then Install Updates.

Figure 3.5: Help file for the print() function.

R Studio might ask you to restart your R session before updating the packages.

3.4 R Markdown

• code chunks
• executing code chunks
• R Markdown reference
• R Markdown cheatsheet

3.5 Some general advice

• naming functions and files
• always use relative links
• so that stuff also works on other people’s computers
• load all the packages at the top of the script
• make sure that a script can be executed from top to bottom
• project management (folder structure)
• don’t write past the vertical rule in code blocks
• learning keyboard shorcuts: Tools > Keyboard Shortcuts Help
• how R handles functions with arguments (order matters if arguments aren’t named)

3.6 R syntax

• base R vs. tidyverse

• two different coding styles

• the pipe

• base R cheatsheet

4.1 Loading packages

The order in which packages in R are loaded matters!

You can refer to functions from specific packages by adding the function name at the beginning. For example, this command would use the select() function from the MASS package MASS::select(), while this command would use the function from the dplyr package dplyr::select().

Always load library("tidyverse") last because it loads a large number of functions that are frequently used.

4.2 Importing data

df.data = read_csv(file = "../../data/top2018songs.csv") %>% 
  mutate(rank = 1:nrow(.))

The quickest way to take a look at your data is to hover your mouse over a variable of a data frame, and press F2.

• importing data cheatsheet

4.3 Data visualiztion

4.3.1 How not to visualize data

We should always take a look at the data first.

include_graphics("../../figures/plots/bad_plot1.png")

Figure 4.1: A not so good plot.

include_graphics("../../figures/plots/bad_plot2.jpg")

Figure 4.2: Another could-be-improved plot.

This second plots reminded me of the following:

include_graphics("../../figures/plots/correlation_aint_causation.png")

Figure 4.3: Correlation is not causation.

Just because two lines look similar, doesn’t mean that anything interesting is going on – it certainly doesn’t mean that the two phenomena represented by the lines are causally connected. For more inspiration check out this site https://www.tylervigen.com/spurious-correlations.

4.3.2 Why you should always visualize your data first

Figure 4.4: The Datasaurus Dozen. While different in appearance, each dataset has the same summary statistics to two decimal places (mean, standard deviation, and Pearson’s correlation).

The data sets in Figure 4.4 all share the same summary statistics. Clearly, the data sets are not the same though.

Tip: Always plot the data first!

Here is the paper from which I took Figure 4.4. It explains how the figures were generated and shows more examples for how summary statistics and some kinds of plots are insufficient to get a good sense for what’s going on in the data.

include_graphics("../../figures/plots/box_violin.gif")

Figure 4.5: Boxplots can be misleading.

ggplot(data = df.data,
       mapping = aes(x = danceability,
                     y = valence)) + 
  geom_point()

ggplot(data = df.data,
       mapping = aes(x = danceability,
                     y = valence)) + 
  geom_point() +
  geom_smooth(method = "lm")

ggplot(data = df.data,
       mapping = aes(x = mode,
                     y = valence)) + 
  stat_summary(fun.data = "mean_cl_boot",
               geom = "linerange") +
  stat_summary(fun.y = "mean",
               geom = "point",
               size = 3)

df.plot = df.data %>% 
  mutate(mode = factor(mode,
                       levels = c(0, 1),
                       labels = c("minor", "major")),
         key = factor(key,
                      levels = 0:11,
                      labels = c("C", "C#", "D", "D#",
                                 "E", "F", "F#", "G",
                                 "G#", "A", "A#", "B")))

ggplot(data = df.plot,
       mapping = aes(x = key,
                     y = energy,
                     group = mode,
                     fill = mode)) + 
  # add individual data points 
  geom_point(mapping = aes(color = mode),
             position = position_jitterdodge(dodge.width = 0.7,
                                             jitter.width = 0.1,
                                             jitter.height = 0),
             alpha = 0.3) + 
  # add the error bars 
  stat_summary(fun.data = "mean_cl_boot",
               geom = "linerange",
               position = position_dodge(width = 0.7),
               size = 0.75) +
  # add the mean data points 
  stat_summary(fun.y = "mean",
               geom = "point",
               shape = 21,
               size = 3,
               position = position_dodge(width = 0.7)) +
  # add the vertical lines
  geom_vline(data = tibble(key = 1:10), 
             xintercept = seq(from = 1.5, to = 11.5, by = 1),
             linetype = 2,
             color = "gray80") + 
  # set title and subtitle of plot 
  labs(title = "Energy for songs with different key and mode",
       subtitle = "Energy represents a perceptual measure of intensity and activity.") + 
  # change the y-axis 
  scale_y_continuous(breaks = seq(0.25, 1, 0.25),
                     labels = seq(0.25, 1, 0.25),
                     limits = c(0.25, 1)) +
  # set the fill color 
  scale_fill_brewer(palette = "Set1") +
  # change the plotting theme
  theme_classic() +
  # adjust the text size
  theme(text = element_text(size = 16),
        plot.subtitle = element_text(size = 12))

4.4 Data manipulation

4.5 Statistics

• linear model lm()
• linear mixed effects models lmer()
• Bayesian models brm() (using library("brms"))

4.6 Saving data

4.7 Help others help you

• making reproducible examples

5.1 Free online books

• R for Data Science: The tidyverse bible.
• YaRrr! The Pirate’s Guide to R: Nice general introduction to R (using mostly base R).
• Learning statistics with R: Introduction to statistics using R (using mostly base R).
• Statistical thinking for the 21st century: Course notes for psychology undergraduate statistics course taught by Russ Poldrack here at Stanford.
• Statistical methods for behavioral and social sciences: Course notes for grad statistica course I teach (notes are not fully self-explanatory though).
• R graphics cookbook: Quick intro to the the most common graphs with ggplot2.
• Data Visualization – A practical introduction (by Kieran Healy): Great resourse for data visualization with ggplot2.
• Fundamentals of Data Visualization: Another excellent resource for data visualization with ggplot2.